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Abstract 

We introduce a random two-matrix model interpolating between a chiral Hermitian (2n + v) X 
(2n+v) matrix and a second Hermitian matrix without symmetries. These are taken from the chiral 
Gaussian Unitary Ensemble (chGUE) and Gaussian Unitary Ensemble (GUE), respectively. In the 
microscopic large- n limit in the vicinity of the chGUE (which we denote by weakly non-chiral limit) 
this theory is in one to one correspondence to the partition function of Wilson chiral perturbation 
theory in the epsilon regime, such as the related two matrix-model previously introduced in [201 121] . 
For a generic number of flavours and rectangular block matrices in the chGUE part we derive an 
eigenvalue representation for the partition function displaying a Pfaffian structure. In the quenched 
case with v = 0, 1 we derive all spectral correlations functions in our model for finite-n, given in 
terms of skew-orthogonal polynomials. The latter are expressed as Gaussian integrals over standard 
Laguerre polynomials. In the weakly non-chiral microscopic limit this yields all corresponding 
quenched eigenvalue correlation functions of the Hermitian Wilson operator. 



1 Introduction 



The application of random matrix theory (RMT) to Quantum Chromodynamics (QCD) first intro- 
duced in [lj has become much more sophisticated in the recent past. Starting from a Gaussian RMT 
with Nf massless flavours [H [2], the so-called chiral Gaussian unitary ensemble (chGUE), several 
milestones have been taken from which we mention only a few, see [3] for recent reviews and more ref- 
erences. After the universality of the microscopic origin limit for non-Gaussian RMT was understood 
[I], the computation of unquenched density correlation functions of the QCD Dirac operator with 
non- vanishing quark masses followed [5j|6], as well as the computation of individual Dirac eigenvalues 
|H II lj. It was understood [10] that the RMT approach is in one to one correspondence to the 
epsilon regime of chiral perturbation theory (chPT) [11] as a limiting theory of low-energy QCD. Cor- 
rections to this regime were computed by including finite volume corrections to the chiral condensate 
£ or to the Pion decay constant F |12j (see [13] and|14j for recent discussions). However, the 
latter only appears in RMT when extending it to a two-matrix model, by coupling to a real [15] or 
imaginary chemical potential |16^ I17j. Otherwise meson correlation functions have to be considered 
in chPT to be sensitive to F, by coupling to the non-zero momentum modes of the Pions. Most 
of the RMT predictions have been checked using Lattice QCD, ranging from checks of the topology 
dependence of algorithms with good chiral properties on small Lattices [18J, up to fully unquenched 
two-flavour simulations leading to realistic values of £ and F [19], to where we refer as well as to [3] 
for more references. 

Very recently a further extension of RMT for QCD has been proposed in [20] in order to include 
also the effect of finite lattice spacing a close to the continuum. The corresponding RMT is again a 
two-matrix model, as it was discussed in more detail in |21j . It is the aim of this paper to investigate 
this RMT further - in fact in a slightly modified version - and to show its integrability by computing 
all density correlation functions. 

The corresponding framework in continuum effective field theory is Wilson chiral perturbation 
theory (WchPT), see [22] for introductory lectures and references. Here, in addition to £ and F three 
further low energy constants W6,7.8 appear to leading order in the infinite volume limit, that have to 
be determined non-perturbatively. The spectral properties of the Wilson Dirac operator were put into 
focus in [23 1, and the epsilon regime in WchPT was first analysed in [24] . The computation of the 
quenched spectral density in the epsilon regime - both for the Hermitian Wilson Dirac operator as well 
as for the real eigenvalues of the non-Hermitian Wilson Dirac operator - was performed in [20J (see also 
[25] for more recent results in the p-regime) . The presentation in [201 [2l] was in fact solely based on 
the WchPT Lagrangian and corresponding generating resolvents, and not on an RMT computation, 
although the two become equivalent in the microscopic limit once the graded or super symmetric 
approach is taken. The extension to the Nf = 1 [26] and to Nf > 1 |27] flavoured spectral densities 
followed also the super symmetric approach. 

Our aim is to partly extend these results and to compute all spectral densities from RMT. The 
advantage of the approach presented here using skew-orthogonal polynomials is that an explicit eigen- 
value basis can be found for the Hermitian Wilson Dirac operator. All spectral correlations then follow 
both for finite- and large- N, where N = 2n + v is the size of the matrices, given the corresponding 
skew-orthogonal polynomials can be constructed. We have succeeded in this program in the quenched 
approximation and for v = 0, 1 as a first step, where v counts the number of exact zero-modes at 
a = 0. We will come back to the open challenges left to our RMT approach in the conclusions. 

Our investigations presented here have a second motivation. So-called parametric density correla- 
tions have been studied in the RMT literature for several decades. The reason of interest is that the 
classical Wigner-Dyson ensembles, the Gaussian unitary, orthogonal and symplectic ensemble (GUE, 
GOE and GSE, respectively) correspond to Hamiltonian systems without (GUE) or with time-reversal 
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symmetry (GOE or GSE, depending on the spin). For that reason the symmetry breaking has been 
studied in terms of a two-matrix model that interpolates between two classes, see [28] for the classical 
papers, or chapter 14 in |29| for a more comprehensive presentation of the GUE-GOE and GUE-GSE 
transition. Several such transitions have been studied since, including the chiral versions of these 
transitions [30] . 

A crucial step is to be able to perform the angular integrations in the term coupling the two 
random matrices. This is done using the Harish-Chandra-Itzykson-Zuber (HCIZ) integral [31J and 
explains, why in all given examples (at least) one of the two classes has to possess unitary symmetry. 
Once an eigenvalue basis is determined the standard (skew)-orthogonal polynomial approach can be 
taken, where in [32j a more general framework for general, non-Gaussian weight functions has been 
developed. 

These two-matrix models also called transitive ensembles include also transitions within the same 
symmetry class, that is from the GUE to GUE or chGUE to chGUE, see e.g. [T7] including their 
extensions with Nf flavours, in the context of QCD with imaginary isospin chemical potential in three 
and four dimensions. As we will show below the RMT corresponding to WchPT in the epsilon regime 
in \20\ [2~T] corresponds (after a minor extension) to the symmetry transition from the chGUE to the 
GUE. Our paper thus serves also to close a gap within RMT by studying this transition class. Because 
of the ubiquity of RMT we expect that applications in totally different areas may follow. 

The presentation of the remaining chapters is organised as follows. In section [2] we present our 
two- matrix model, show its relation to WchPT in the epsilon regime (see also appendix lAl for details) 
and discuss the relation to the RMT proposed in (20JET]. In section [3] we derive the joint probability 
density function for the Hermitian Wilson Dirac operator for a general Ay-flavour content and general 
v > 0. The solution of our two-matrix model for all density correlation functions is presented in section 
[3] in terms of skew-orthogonal polynomials for Nf = and u = 0,1, see also appendix ICl for further 
details. These are given in terms of Gaussian integral transforms of standard Laguerre polynomials 
as they appear in the chGUE [lj. The microscopic large- N limit at the origin close to the chGUE is 
then presented in section [5] illustrated with several examples, before we present our conclusions and 
open problems in section El 



In this section we first introduce our matrix model and explain why it describes the transition between 
the different symmetry classes of random matrices from the chGUE to the GUE. Then we point out 
the relation to the related Wilson chiral RMT previously introduced in [20J as well as to WchPT in 
the epsilon regime. 

We want to compute the eigenvalue density correlation functions of the following Hermitian Wilson 
Dirac operator: 



Here W is a rectangular matrix of size n x (n + v) with complex elements, without further symmetry 
restrictions. The non-negative parameter v > is related to the number of zero-eigenvalues in the 
limit to the chGUE. The real parameter m denotes the quark mass, as will become clearer below. 
This first part of 05 is what we are used to in the chGUE (up to the sign in m). The second part 
is a quadratic Hermitian matrix of size N = 2n + v with complex elements. It corresponds 
to the GUE part of breaking chiral symmetry. 



2 The Random Matrix Model 




+ H . 



(2.1) 
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The matrices W and H are distributed with the following Gaussian probability measures: 



Pi(W) 
P2{H) 



exp 
exp 



2(1 
1 
4^2 



-TrWW* 



a G [0, 1] 



(2.2) 



We expect that the simplest, Gaussian choice that we have made here is not important in the large-n 
limit due to universality. We can now define the following Nf flavour partition function 



zi Nf \m;{z f };a) 



dHdW Y[ det[0 5 + z f l N ]Pi(W)P 2 (H) , 
/=i 



(2.3) 



where we have inserted a product of Nf characteristic polynomials of real arguments Zf, in addition 
to the weight functions. The integration is over all independent matrix elements of W and H. 

The role of the real parameter < a < 1 is to interpolate between two different symmetry classes. 
For Nf = and m = the limit a — > leads to a delta function in all matrix elements of H , reducing 
eq. (j2.3j) to the chGUE. In the opposite limit a — > 1 we obtain a delta function in all matrix elements 
of W, and we are lead to the GUE of size N. For m / we interpolate between the eigenvalues 
of the Dirac operator with finite quark mass m and the GUE coupled to a fixed external field with 
eigenvalues ±m, as it was considered for example in [33] ■ 

In |20j a very similar Wilson chiral RMT was introduced which we will label by II. Starting from 
a non-Hermitian Wilson Dirac operator D\v, 



D w 



aA W 
-W* aB 



(2.4) 



with W as above and A = n x n and B = B< (n + v) x (n + v) Hermitian matrices, respectively, 
the Hermitian Wilson Dirac operator D§ in [20] was obtained by multiplying Dyy + m with 75: 



D R 



1, 







(D w + mliv) 



ml r 



W 
-mlr, 



+ 



aA 






aB 



-i n+u j v iyt 

which is Hermitian. The corresponding partition function with Gaussian weights reads 



(2.5) 



z S } ( m ;{ z /}' a ) ~ 



Nt 



dHdW det[D 5 + z f t N ] exp 
/=l 



1 



Tr(yl 2 + B 2 + 2WW r ) 



(2.6) 



Because of symmetry the minus sign in the lower right corner in B in eq. (|2.5p can be absorbed by 
shifting B — > —B in the integrand. 

Both matrix models Zv and Z^f, enjoy the same relationship to WchPT in the epsilon regime 

z (see appendix [A] for the non- 

(2.7) 



in the following large- iV~ limit, stated here for equal arguments Zf 
degenerate case): 



lim zi Ni \m; z;a) 



lim Z^\m;z;a)=Z^\m,z,d 

N—too ' 



where 



i Nf \m,z,a) = ( dUdeb[U]"e* 

Ju(Nf) 



2 P svTr(c/+(7t)+i2pSvTr(c/-(7t)-ap 2 vvi/8Tr((7 2 +c/t2) 



(2.8) 
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Here we have to match the large-iV with the large- U infinite volume lining] 

m = m p YAf = mVN , I = z p Y,V = zVN , and a 2 = a 2 p VW s = ^a 2 N . (2.9) 



While the equivalence eq. (|2.7p for the second RMT eq. (|2.6p was shown in [20] we will show this 
equivalence for our model eq. (|2.3I) in appendix lAl 

Eq. (|2.8p . the partition function of WchPT after Fourier transform, usually contains two further 
low-energy constants Wq and Wj, apart from the infinite volume chiral condensate S, and W% which 
encodes the effects from a finite-lattice spacing to order 0(a 2 ), which we displayed. Because Wq and 
W-j can be switched on at the expense of one additional Gaussian integral each |21j . we will only 
consider Wq = W7 = in the following. The parameter m p denotes the standard (equal) quark-mass 
term coupling to vpip in field theory, whereas the z p denotes an additional source terms coupling to 
^75^ which will be convenient later. 

Because of the matching eq. (|2.7h we will assume that both RMT eqs. (12, 3h and (|2.6h as well as 
WchPT eq. (|2.8I) are in the same universality class, also regarding all eigenvalue density correlation 
functions. 

The matching between the two matrix models can be made more explicit on the level of matrix 
elements. By redefining H/a — > H and Wj \J\ — a 2 — > W in eq. (|2.ip we obtain to leading order 

m ( mt n W \ ( aA att + ±a 2 W\ , 



where 



The advantage of allowing for extra, off-diagonal matrix elements in matrix H in our model is that 
we can do a proper change of variables from {W, H} to {W, P5}, keeping the same number of degrees 
of freedom. This allows us to go to an eigenvalue basis as we will show in the next section. 

Furthermore, the choice of parametrisation in a in our model is more convenient when studying 
the symmetry transition from the chGUE at a = to the GUE at a = 1. Compared to that in D§ eq. 
(|2.5p the transition starting at a = corresponds to the chGUE too, then leading via a = 0(1) of a 
full N x N GUE to a decoupling for a>l into two GUEs of sizes n and n + v respectively. 



3 The eigenvalue picture 

In this section we will derive an eigenvalue representation of our partition function eq. ()2.3p . It can 
be obtained by the following change of variables from H to 05, 

-JL)- <"> 

then diagonalising 05 and W, using the HCIZ formula, and then finally integrating out the W- 
eigenvalues. 

1 In |20l 121] all matrix elements are rescaled by y/N to have compact support at TV = 00. This will modify these 
scaling relations on the RMT side by \/~N for z and m. 
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The Jacobian of the change of variables eq. (13.ip is a constant and we get for the partition function 



zl Nf \m;{z f };a) ~ / dfcdW J] det[0 5 + z f l N }P 1 (W)P 2 [$, 

J f=i V 



?7ll r 



w 



-ml n j 



e 4o-* 



AT, 



x / ciWexp 



<i05 det[0 5 + z/l/v] exp 
1 



4a 2 



Tr[0 



2a 2 (l-a 2 



-Tr 



-I Tr 

2a 2 



^ , ml n W 
P5[ -mln+u 



(3.2) 



Now we can diagonalise Ij)^ by a unitary transformation U £ U(N), 

]p 5 = UDU ] , D = diag(di, . . . ,d n ) , dj G R , (3.3) 
and the chiral matrix W by a singular value decomposition (or equivalently a diagonalisation of WW^): 



W 



uYv 44> 



mi n W 
-mtn 



ut 
it 



y 

-mt n+v 





u 



(3.4) 



where Y is a rectangular n x (rt + z/) matrix with real positive diagonal elements {yi, . . . , yj > 0. 
Note the different numbers N = 2n + v and n of these eigenvalues dj and , respectively. Including 
the Jacobians of these transformations, which contain the standard Vandermonde determinants, 



N 



Aiv({4) = - ^ > A «(^ 2 » = iw - y" 

j>l j>l 



(3.5) 



we get 
zl Nf (m;{zf};a) 



e 



TV N f 



Udd.Hid.+z^e-^ANdd}) 2 / \[dy b yl v+l 



'i=i /=i 



xA n ({?/ 2 }) 2 / dUdudv exp 



6=1 



2a? 



-ml. 



where we have used the invariance of the Haar measure under the transformation 

ft 



U 



u . 



u 

The unitary integrations can now be carried out and we obtain from the HCIZ integral 



(3.6) 



(3.7) 



dU exp 



1 

2a7 



Tr 



UDU ] ( 



ml. 



V Y 



Y 

-mL 



det 



l<i,i<iV 



(3.8) 



(mi. y 
F J -ml n+p 

relate to the eigenvalues y& in order to partly cancel the Vandermonde determinants. From the 
solutions of the characteristic equation 



A N ({d})A N ({w}) 

which we will now 



det 



Al 



N 



ml n Y 
Y T -ml n+u 



(\ + mYl[(x 2 -m 2 -y 2 b 



N 



6=1 



H (A - Wj ) ,(3.9) 
i=i 
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the eigenvalues Wk are given bjH 



m 2 + y 2 , j = 1, 



,n 



w 



m 2 + y 2 , j = l,...,n 



Wj + 2n 



-m + ej ,j = l, 



, v ■ 



(3.10) 



Here we have lifted the degeneracy of the last v eigenvalues by adding small parameters e% , . . . , e u 
which will be sent to zero later (for v = 1 this is not necessary, but we will keep e\ ^ for symmetry 
reasons). Making the e-dependence explicit in the Vandermonde we can write 



In v 



A N (w) = A u ({e - m})Y[Y\( e h ~ m ~ Wl) Y\_( w i+n ~ W j+n )(Wi - Wj) Y[ ( w o+n 

1=1 h=l i>j o,g=l 



(-^H^nil^-^-rf) A n ({y 2 }) 2 Jl Jm^+yf , (3.11 



1=1 h=l 



i=l 



after some algebra. The first Vandermonde determinant A„({e}) will cancel the zeros from the numer- 
ator in eq. (|3.8[) when we take the limit of degenerate eigenvalues — > 0, whereas in the remaining 
factors in eq. (13. lip this limit is smooth. 

Using induction one can easily obtain the following limiting result. For v = 3 we have to Taylor 
expand the next column to one order higher, when taking £3 — > e±, etc., and so we end up with the 
following result to leading order 



1 



e k ^o A„({e}) 



u-l 

n 

j=0 



3 ,xiW2n „xi(—m+ei) 



e x 1 W2n e —xim x\e~ Xim 



^,x N wi 



JCNW2n Q — X N m XNe~ XNm 



Jf-1, 



™V—l -x N m, 



(3.12) 



where we have defined X{ = ^jdi for i = 1, . . . , N for later use. 

We can now write the full answer of our partition function eq. f|3.6|) in term of the two sets of 
eigenvalues only: 



Zu f \m;{zf};a) ~ e 



^ 2 (2n + JV(l-a 2 )) rOO N 



Aa 2 (l-a 2 ) 



Nf 



II dd i YiS d i + Z f^ e ^ A N({d}) / \~[ du b e 



' i=i /=l 



. . . e 



1 2a z (\-a 2 ) 



6=1 



x\ e 



...Xi 



1XNU1 p x N u n 



-xnu\ 



. . . e 



„u— 1 — xjym 
• • • -kjv c 



(3.13) 



after dropping m-independent constants and changing variables Uj = x jm 2 + y 2 . For v = the last v 
columns that are independent of Uj have to be dropped. 



2 From this it can be seen that at a = there is an accumulation of eigenvalues at — m. Note that our convention 
differs from |20I|2T1. 
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We will now apply a generalisation of the de Bruijn integration formula in order to integrate out 
the set of variables Uk- It is given by [33] (see appendix C.2 therein) 



II / du i 

q=l J 



det 

l<b<2n+v; l<j<n; l<i<v 



{<t>b{Uj),^b{Uj)} A U 



(") 



u(u-l)/2 . 



ni Pf 

Kb,c<2n+u; \<%<v 



{J duw(u)[4> h (u)il) c (u) - 4>c(u)^b(u)]} A bi 
-At 



(3.14) 



Here A^ = d l b 1 e d b m / 2a2 \ s a matrix of size N x v. When v = it is absent and we are back to the 
standard de Bruijn formula. In our case we have 



zi Nf (m; {zj}; a) 



exp 



xPf 



Nm 2 



2a 2 (l - a 2 ) 



oo 2n+i> 

n n (4? + z f) ex p 



4a 2 



A 2n+ „({d}) (3.15) 



l<*J<2n+;/; l<g<^ 



where we have defined the antisymmetric weight 



Fix) 



e 8a" 



erf 



8a 2 

oo 



+ 



+ erf 



x 



VT 



/?7 



3a 2 ^2a 2 (l - a 2 ) 



(3.16) 



cfot e 2o2(l- a 2) 



xu 



2a 2 



We have also used the antisymmetry of erf(x) = Jq dse s2 = — erf(— x). Eq. (13. 151) is the partition 
function in terms of the joint probability distribution functions (jpdf ) of its eigenvalues. 

Let us give the two examples that we will solve explicitly in the following sections. For Nf = v = 
we have 



30 2n 

Y[ ddj exp 

00 3=1 



4a 2 



A 2n ({d})Ph< iJ <2n[F(d j -di)] . 



(3.17) 



Here we have omitted all constants that will drop out in expectation values and density correlation 
functions to be considered in the next section. For v = 1 and Nf = we have 



i£f°W) = 



oo 2n + l 

Y[ ddj exp 



L 4a 2 



A 2n+ i({d})Pf 



l<i,j<2n+l 



F(dj — di) e~ 



_ d j n 
-e ^ 







. (3.18) 



Let us add one remark on the universality of our model. In order to be able to use the HCIZ 
formula it was crucial that we started with a Gaussian distribution of matrix elements. However, after 
having arrived at an eigenvalue representation eq. (|3.15p we could take this as a starting point and 
we would obtain the same results below, if we were to replace the Gaussian distribution exp[— d 2 /4a 2 ] 
by a more general weight function w(dj). The same general framework to compute density correlation 
functions to be presented in the next section could be applied. 
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4 All density correlation functions for finite- TV 



The partition function eq. (|3.17p is very much reminiscent of the one for the GUE-GOE transition in 
[28] . apart from the different function F(x) inside the Pfaffian. We can therefore apply the general 
method of solving this two-matrix model developed in [32] for general weight functions, which applies 
equally to eq. ()3.18|) with v = 1. From now on we will restrict ourselves to v = 0, 1 and Nf = 0. 
The /c-point density correlation functions that we will determine are defined as follows for u = 

N 



Pk(dx, ... ,dk) 



1 



/OO J ' 

dd k +i ■ ■ ■ dd N JJ w(dj) A N ({d})Ph< i;j < N [F(dj - di)] , (4.1) 



with even N = 2n, and for v = 1 as 
1 N\ 



Pk l (di,--- ,d k ) 



N 



dd k+ i ■ ■ • dd N ^[w(dj)&N{{d})¥ii<ij<N 

3=1 



F(dj-di) f{di) 



-/to: 



with odd N = 2n + 1 . Here we have used the abbreviations 

w(x) = exp[— x 2 /(4a 2 )] , f(x) = exp[— mx/(2a 2 )] . 



(4.2) 



(4.3) 



The expressions for the solution in terms of three kernels will hold for a general weight functions 
w(x), f(x) and F{x) though. Here and in the following we will drop the label v = and only index 
v = 1. Furthermore we also suppress the dependence on the parameters m,a in the arguments for 
simplicity. 

The correlation functions defined above can be expressed as the Pfaffian of the following matrix 

[32] 

In(di, dj ) Sn(di, dj ) 
-SZ(dj,di) -D-idudj) 



Pk(di,. . . ,d k ) = Ph<i,j 



<k 



(4.4) 



given in terms of three kernels. 

Let us begin with the case v = where we have [32] 



S n (x,y) 



j2^1U 2j _ 2 ( x )R 2j -i{y) - faj-i(x)R2j-2(yj) , 

■ i 7 — 1 

^ wxjwjy fji 2 . 2 ^ji 2jl ^ _ R 2j _ l (x)R 2 j-2(y)) , 

n 

-^2 (hj^^^j-iiy) - <j)2j-i(x)(f>2j-2(y)) - F(x-y) . (4.5) 

• 1 r i— l 



In our specific case of weights eqs. (13.160 and ()4.3p these kernels are given by the following monic 
skew-orthogonal polynomials Rj(x) = x 3 + . . . and their integral transforms, 



R 2 j(x) 
R 2 j+i(x) 

4>j(x) 



j! 2J(1 -a 2 ) 
(-l)iV5F 



2\j /-oo 



ds e s La 



[x + 2iasY 



m 



d r 



-2a 2 u;(x) _1 — • 
ax 

i!2J'(l-a 2 y /- 00 



(-l)i0F 

OO 

w(y)F(x - y)Rj{y) 



w(x)R2j(x 

2 

cis e~ s (x + 2ias)Lj 



2(1 - a 2 ) 

x + 2ms) 2 



(4.6) 



m 



2(1 -a 2 



(4.7) 
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Here we have used an integration by parts to simplify the expression for the odd polynomials. The 
polynomials have parity, Rj(—x) = (— as can be easily seen from the integral representations. 
These skew-orthogonal polynomials satisfy the following relations as it is shown in appendix \C\ in 
detail: 

(i?2j) R2l+l) = r j$jl ! 
(R 2 j, R21) = = (i?2j+i, R21+1) (4-8) 

with the skew product defined as 

(h,g)= / dxdyw(x)w(y)F(x - y)h(y)g(x) . (4.9) 



The squared norms read 

Tj = (R 2j ,R2j+i) = 8V27ro2(l - a 2 ) 2j+2 2 2 i(j\) 2 e , (4.10) 

are also derived in appendix The quenched partition function can now be expressed in term of 
these norms and is given by 

n-l 

Z<%*>=mUr j . (4.11) 

j=0 

Because of the special form of our polynomials the kernels can be simplified. We start with the 
kernel D n (x,y), the other two easily follow from the relations 

/CO 
dz F(x - z)D n (z,y) , 

/DO 
dz F(x- z)S n (x,z) . (4.12) 
-00 

Inserting eq. (j4.6|) into the definition we obtain 

m 2 

x 2 +y 2 e 2(l-a*) j-OO _ 2_ 2 

D n (x,y) = e / dsdre s r \(y + 2ias) — (x + 2iar)] 

V ' 8vrV2vra 2 (l - a 2 ) 7-oc 



^ / (x + 2iar) 2 -m 2 \ / (y + 2ms) 2 - m 2 \ 

^ Lj 1 2(T^j J Lj V 2(1^ ) • (413) 



x 

3=0 



Using the Christoffel-Darboux identity 

n fL j{ X)L i{ Y) = -„ ^(X)L^ l{ Y)-L (Y)L n ^X) ^ 

3=0 

we obtain the following final expression 



iV eW^nV] - a 2 f°° , , _ c 2_^2 



D n (x,y) = e~~^ j= / dsdre~ s - r (4.15) 



— 00 



x 



j ( (y+2ias) 2 —m?\ j ( (x+2iar) 2 — m 2 \ j ( (y+2ias) 2 —m 2 \ 
n \ 2(l-a 2 ) J ^n-l y 2(l-a 2 ) ) ~ ^-1 \ 2(l-a 2 ) ) ^ n \ 2(l-a 2 ) ) 



(x+2iar) 2 — 



(y + 2ias) + (x + 2iar) 
10 



The simplest examples for the correlation functions at finite- iV are the spectral density and the spectral 
two-point function. These are given by 



pi(x) = S n (x,x)= / dz F(x - z)D n (z,x) , 

J — oo 

P2{x,y) = S n (x,x)S n (y,y) + I n (x,y)D n (x,y) - S n {x,y)S n (y, x) , (4.16) 

with the insertion of eq. (|4.15|) and (|4.12|) . see eq. (|4.25p below. This ends the general solution of our 
model for finite- N for v = 0. 

For v = 1 the skew-orthogonal polynomials can be expressed in terms of those for v = 0. Following 
|32j we introduce new polynomials 

m =1 (x) = Rj{x) - -JI-Rn^x), j = 0,l,'~ ,N-2 , (4.17) 

and for the largest polynomial 

R v £-x{x) = Rn-i(x) , (4.18) 

with coefficients 

/oo 

dx w(x)f(x)Rj(x). (4.19) 
-oo 

Moreover we define the integral transforms of the polynomials as in the previous case eq. (|4.7p 

/oo 
dy F(x-y)w(y)R» =1 (y). (4.20) 
-oo 

The kernels for v = 1 to be inserted in eq. (14. 4p read |32| 

w{y) 



S v = l (x, y) = — (^W^H^ " ^(*)HgtM) 
1 



+ f(x)w(y)R» N =\(y), 

SN-l 

w(x)w(y) 



Dt\x, y) = J2 (R» 2 =} 2 (x)R» 2 =_\(y) - RZTMx^My)) , 

n 
3=1 

+— ttW-lWttv) - <t>N=My)f(x)) - Fix - y). (4.21) 

SN-l 



The partition function is also evaluated as 



n-1 



4=r 0) = ^^-in^-- ( 4 - 22 ) 

j=0 

As a last step we still have to determine the new coefficients Sj. We can readily find 
s 2j = / dx e-< x+2m V^R 2j (x) 



HP loo 3 { 2(1 -a*) 

(_l)i j! 2 J +1 (1 - a 2 )^ a e m2/(4a2) , (4.23) 
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for the even coefficients, using an argument as in eq. (|C.19j) . and 

/OO A r 

-oo 

/"OO 



oo 



= -m s 2 j , (4.24) 

for the odd coefficients. This ends our general solution for v = 1. 

In order to check that the limit a — > correctly reproduces the chGUE let us spell out the densities 
explicitly, where we start with v = 0. Here we use eq. (|4.13p rather than eq. (|4,15p and we obtain 



px(x) = dz 



oo 



g 8a-! r 8 ~ r 2(l-a 2 ) 

8^^271(1 - a 2 ) 



,. (x — z)(l — a 2 ) + 2m\ / (x — z)(l — a 2 ) — 2m\ 
erf = + erf 



2av/2(l - a 2 ) / V 2(1^/2(1 - a 2 ) 



/•°°,, _ s 2_ r 2 ^ /(z + 2mr) 2 -m 2 \ /(x + 2ias) 2 - m 2 \ 

(4.25) 



Because of the delta-function that we obtain from 



1 (z + x) 2 

lim — e ~~^r~ = S(z + x) , (4.26) 



a^o 2aV2vr 

we get the following from the error-functions (cf. IB.9|) 

J 2x(l - a 2 ) ±2m\ . , A . 

hm erf — - — ; = sign(x ± m) . (4.27) 

^ 2a^2(l-a 2 ) ^ V ' K ' 

The resulting expression we can rewrite as 

x(sign(x + m) + sign(x — m)) = |x|2©(|x| — m) . (4.28) 

In the second line in eq. ()4.25p the integrals decouple and we obtain a simple sum of Laguerre 
polynomials. The final answer is thus 

2 2 n ~ 1 / T 2 m 2\ 2 



i=o 



limpi(x) = |x|e x a + ™ ( ) 0(|x|-m) (4.29) 



W pf GC/i? (Vz 2 -m 2 ) 6(|x| - m) , (4.30) 
yx 2 — m 2 ^ ' 



which corresponds to the shifted density of the chGUE for finite-n, due to the extra mass in our Dirac 
operator definition eq. (|2.ip . 



ra-1 

,,2 



pi hGUE (y) = \y\e-y 2 / 2 Y,L l (y 2 /2) 2 . (4.31) 



z=o 



As an illustration we have plotted the densities in figure [TJ For comparison we also display the 
finite-n density of the chGUE eq. (|4.3ip in the limit a — > 0. Switching on a > the smoothening of 
the hard wall provided by the Heaviside-Theta functions is nicely seen in figure [TJ 
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Figure 1: Example for the v = quenched density at finite n = 4 at different masses: m = 0.4 (top 
left plots) and m = 1 (top right plots): the top black curves correspond to the shifted chGUE density 
eq. (|4.29p at a = (which is clipped for m = 1) ending at a sharp cutoff for m > 0, see eq. (|4.3Qj) . 
The density pi(x) eq. (|4.25j) is plotted at various values of the deformation parameter a: at a = 0.1 
(top smooth blue curve), a = 0.3 (middle smooth red curve), and at a = 0.9 (lowest yellow curve) 
which is already indistinguishable from the GUE curve. Because the spectral density is symmetric 
due to v = we only show the positive eigenvalues. The lower plots show the same curves at m = 0, 
the chGUE-GUE transition. Here a = and a = 0.1 are already indistinguishable. 

Let us now spell out the density for v = 1 explicitly and perform the limit a — > there. In view of 
the relation between the polynomials eq. (|4.17p it is useful to express the kernels for v = 1 in terms 
of those with v = as well. We begin with the simplest kernel containing only polynomials and no 
integral transforms. After some algebra we obtain 

W -1/ x / \ rr? i r roo / (y + 2ias) 2 - m 2 



2\ 



8aTT^2iv{l - a 2 ) U-oo V 2(1 - a 

x J°° dr e- r \x + liar + m)X^ x — ) ~ (z^y)} , (4.32) 

where we used the following identity for generalised Laguerre polynomials [35] 

Y J L j {z) = L^ 1 {z) . (4.33) 

For the kernel determining the density we obtain after inserting the definitions 

■ 2±, -- L " a 1 r °° / (y + 2ias) 2 — m 2 \ 

(4.34) 



QU=1 



(x,y) = / dzF(x- z)D u n - 1 (z,y)+e 53 — / ds e~ s L 
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These two equations together with eq. (I4.15P then determine the spectral density for v = 1 shown in 
figure [2j 

P T\x) = S v =\x,x) . (4.35) 




Figure 2: Example for the v = 1 quenched density at finite n = 4 at different masses: m = (left 
plots) and m = 1 (right plots). Because for m = the density is still symmetric we only plot the 
positive eigenvalues. The black curves correspond to the (shifted) chGUE density at a = 0. It can be 
nicely seen how the delta- function at x = —m gets gradually smoothed out, where we show a = 0.1 
(blue curve) and a = 0.3 (red curve). 



We can now perform the limit a — > as a check. The last term in eq. (I4.34D gives a delta-function, 
8(x + m), the exact zero-eigenvalue of the v = 1 chGUE shifted to x = —m. Using the result eq. 
(|4.29h and the above results for the limit at v = we obtaiiil 



limpi 1 (x) 

o->0 



\x\ e 2 



n-1 
3=0 



2 2\ 2 

x — m x 



G(|x| — m) 



\x\e 2 + 2 L n 



i i 
x — m 



L 



(i) 

n—1 



2 2 

X — m 



2 J n ~ l \ 2 

This can be seen as follows to be equivalent to the density of the chGUE at v = 1 

n-1 

Pl,v=l 



Q(\x\ -m) + 8{x + m). (4.36) 



.chGUE iy) = ^3 e -yV2- 



1=0 y ' 



shifted as in eq. (|4.3Q|) . When using further identities [35] 

XLf ) (X) = (l + l)(L l (X)-L l+l {X)) , Lf ) (X) = L l (X)+L { l l \(X) 



(and setting L_i = here), we obtain the desired result, with X = \(x 2 — m 2 ), 

n—1 ^ n—1 

x Ettt^ (1) w 2 = E L / (1) w(^w-^ + iw) 

1=0 1=0 



1=0 

n-1 



£ (MX) + L« - ^L z (1) (X)L m (X) 



n-1 



Z=0 

n-1 



«=0 



^LKX) 2 - L^Ux)L n (X) . 



1=0 



(4.37) 



(4.38) 



(4.39) 



3 From this we might speculate that the second term in eq. (|4.34[) leading to 5(x + m) describes the distribution of a 
single eigenvalue at finite- TV. 
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5 The microscopic weakly non-chiral large-N limit 



The microscopic large- ./V limit is defined by the following rescaling of variables 

( S VE., * S \-J2na. (5.1) 



These variables are then to be identified with the variables in WchPT as given in eq. (|2.9p . Here 
we send the eigenvalues to zero and n — > oo such that x remains finite, hence we are in the vicinity 
of the origin. The masses m and z are rescaled likewise. Also the parameter a corresponding to the 
effect of finite-lattice spacing in WchPT is sent to zero and n — > oo such that a remains finite. Thus 
we are in the vicinity of the chGUE ensemble which has chiral symmetry. For this reason we call this 
microscopic origin limit weakly non-chiral. 

We begin with the asymptotic formulas of the partition function. Because we know the (Nf = 1)- 
flavour partition function for finite- N and general v > we compute its asymptotic first. It reduces 
to the asymptotic for the even polynomials R 2 j at v = 0. Based on the standard limit |35j 

lim n-"LM(-) =x-^ 2 J u {2^) , (5.2) 

n— >oo \nJ 

we obtain the following expression: 

lim ^(det[x + 5 ]) = lim (1 - a 2 y [°° *L e~* \x + 2ias - mfL^ ( {x + ^ " ^ ) 
j^oo 23 j\ n^oo J_ OQ ^/vr J \ 2(1 - a z ) J 

r ^ f ^r^ + f a) Vi,(v ^-(x + fea) 2 ) ). 

J-oo V 71 \ m + x + Aisa J \ J 

(5.3) 

First we have divided out part of the normalisation that cancels with the norms rj in the sums inside 
the kernels. Second we have given the result for j = tn with t £ [0, 1] as we will need these later when 
replacing the sum inside the kernel by an integral. Last we have used that due to the rescaling the 
prefactor turns into an exponential, 

(1 - ja 2 /j) j -> exp[-ja 2 ] = exp[-2ta 2 ] . (5.4) 

The result eq. (|5,3p matches that of the partition function eq. (|A.7p as it should. 

The limiting even polynomial R2j{x) for u = is obtained by simply setting v = everywhere in 
eq. (|5.3p (due to parity there is no need to switch x — > —x under the determinant): 

lim \-Lr 2j ( x ) = e~ 2ta / ^e~ s U^t(m? - (x + Aisaf) 
nj->oo 23 j\ J_ OQ Vvr V 



n,j—¥oo 

-2ta 2 



(5.5) 



The computation of the limiting odd polynomials R2j+i(x) for v = is now a simple consequence. 
Starting from the explicit form eq. (|4.6p we obtain 



lim V2^^-R 2j+1 (x) = e -2^ 2 r *L e -s 2 (£ + 4i S a) Io ( ^(^2 _ p + 4^)) , (5.6) 

nj^oo 23 j\ J_ OQ yJlT V / 



after an integration by parts. As a check the asymptotic of the odd polynomial is an odd function in 
x. Note that both prefactors from eq. (|5.4p will cancel with the limiting norms rj. 
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Figure 3: The microscopic density eq. (|5.8p for rh = 1 (top left plots) and rh = 3 (top right plots) 
at various values of a = 0.1 (middle blue curves) and a = 0.25 (lower red curves). For comparison 
we also show the corresponding chGUE density (d = 0) eq. (|5.1ip shifted as in eq. (|4.30p (top black 
curves). For the chGUE-GUE transition at rh = (bottom plots) the difference to a = 0.1 or a = 0.25 
is almost indistinguishable, apart from the origin. 



For the limiting weight F(x) from eq. (13.16P we get a Gaussian times the error-functions: 



lim F(x) = exp 



x 2 



32a 2 



. x m 
erf ( — ^= + 



4V2S2 2VW 



+ erf 



m 



4V2S2 2v / 2d 2 



(5.7) 



Collecting all results we obtain the quenched v = microscopic density ps(x) = limn^oo -^=p 1 (x): 



oo 



e 32a 

dy 



- x 32V2a 2 vr 
l 



dsdr _„2_„2 
e 

7T 



4V2a 2 y v 4V2a 2 
x / dt(x — y + 4id(r — s))Iq (y/ t(m 2 — (y + 4isd) 2 )^ Jo (\/ t(m 2 — (x + 4ird) 2 )J 



e 32a2 



16\/2a 2 7r 



.,rf(^±J^Uerff^^) 



4V2a 2 



4V2a 2 



dsdr _„2_„2 
e 

7T 



m 2 — (x + 4ird) 2 Io ( yxrh 2 — (y + 4isd) 2 I Ii I a/ to 2 — (x + Air a) 2 \ — (x y) 



x + y + 4m(s + r) 



(5.8) 



It is shown in figure [3l Here we give both forms, before and after using the Christoffel-Darboux 
identity which has its large- iV correspondence in 



f l , r, , r, 2XUY)h{X) -2YUX)h 

J dtI (XVt)I (YVi) = ° l ' X 2_ Y 2 



00 



(5.9) 
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Unfortunately we have been unable to check analytically that eq. (15, 8h is equivalent to the corre- 
sponding density given in [20]. There it is given as the discontinuity of the imaginary part of its 
resolvent. 

As a check we can again take the limit a — > 0, and we obtain from the first form in eq. (|5.8p 



lim ps{x) = ^-©{\x\ — rh) [ dtJ® ( y / t(x 2 — rh 2 )) 
a^O 2 Jq \ J 



Q(\x\ — rh) Jq ( \/ x 2 — rh 2 ) + J\ 



x 2 — rh 2 



It equals the shifted microscopic density of the chGUE at v = after the shift in eq. (|4.30p 



„chGUE ( 

Ps 



(x) 



J (x) 2 + Ji(x) 



(5.10) 



(5.11) 



We now turn to the limiting density for v = 1. Starting from eqs. (|4.32|) and (|4.34|) it is useful to 
write it as the limiting density at v = which we have already determined, plus a correction term. 
Collecting all terms and using the asymptotic expressions from above we arrive at 



Ps~ (£) = Ps(x) + 



(x + zf 
e 32a' 2 

dz — - — -= 

oo 16a-7r\/2-7r 



(iH x-z + 2m ]+e ^x-z-2m^ 



ds e 



+ Aisa) 2 



dr e 



rh + z + Aiar 
\l rh 2 — (z + Aiar) 2 



x h(y 

l (x + m) 2 

+ - — re lea 2 



rh 2 



(z + Aird)'- 



(X <R- Z) 



Aird 



which we illustrate in figure HI 



dse s2 Io[\/rh 2 — (x + Aisd) 2 ), 



(5.12) 




Figure 4: The microscopic density eq. (15.120 for rh = (left plots) and rh = 3 (right plots) at various 
values of d = 0.3 (red curves), a = 0.1 (blue curves), and d = (black curves) for comparison. For 
the chGUE-GUE transition at rh = the density is still symmetric and we only display x > 0. For 
rh > the delta-function is shifted to the left, and we can nicely see its increasing broadening from 
d = 0.1 to a = 0.25. 



Once more we can perform as a check the limit d — > 0, and we obtain 

x 



lim p v s (x 

a— >0 



lim ps{x) 



a-5>0 



X* — 777/ 



Jo [v x — rh 2 J\ 



x 2 — rh 2 



^ 0(|x| — rh) + 5(£ + rh) . 

(5.13) 
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The matching with the density of the chGUE at v = 1, 

pf G v Ef{x) = ^(Ji(x) 2 - J (x)J 2 (xj) , (5.14) 
shifted as in eq. ()4.30p can be seen as follows: 

pt G ™(x)-pt GUE (x) = -^j (x)(j 2 (x) + j (x)) 

= -JopOJxpO . (5.15) 

Here we used a Bessel identity and the abbreviation X = \J x 1 — m 2 . The delta-function in eq. (|5.13p 

is just the exact zero-eigenvalue at x = —rh, originating from the last term in eq. (|5.12p which is 
non- Gaussian. 



6 Conclusions 

In this paper we have introduced a parameter dependent Gaussian random two-matrix model that 
interpolates between the chGUE (a = 0) and the GUE (a = 1) symmetry class (possibly shifted by 
a constant matrix). At the same time this model describes the spectral properties of the Hermitian 
Wilson Dirac operator when we take the microscopic weakly non-chiral limit at the origin. Here 
the rescaled parameter a ~ a^/n is interpreted as the effect of finite lattice spacing. It is properly 
incorporated in a continuum effective field theory as Wilson chiral Perturbation theory in the epsilon 
regime to order 0(a 2 ). 

We have completely determined all spectral /c-point density correlation functions when starting 
from a chGUE of size (2n + v) with v = or u = 1 exact zero modes, both for finite n and in the 
microscopic large-n origin limit at weak non-chirality. They are given in terms of the Pfaffian of a 
2x2 matrix kernel, containing three individual kernels as building blocks. Each kernel is expressed 
as a sum over skew-orthogonal polynomials (and their integral transforms) depending on v = or 
1, which we explicitly constructed in terms of Gaussian integrals over Laguerre polynomials. In the 
large-n limit these 3 respective kernels contain in their final form a fixed number of 2 up to 4 Gaussian 
integrals over modified Bessel functions, respectively. The Pfaffian of the matrix kernel is thus easily 
evaluated for any fe-point function. 

Let us compare to what is known in the literature. In [20|, [2Tj the super symmetric method was used 
to compute the microscopic spectral density (k = 1) directly from the chiral Lagrangian for arbitrary 
u, the so-called index. In |26| this was extended to include one (Nf = 1) extra quark determinant. 
This method rapidly becomes cumbersome when adding more flavours Nf > 1 or increasing k, as the 
number of integrations increases with the dimension of the auxiliary supermatrix from the Hubbard- 
Stratonovich transformation. This problem was partly circumvented in [27] by using the graded 
super eigenvalue method, where explicit examples where provided for the density k = 1 up to Nf = 2, 
and for the two-point density k = 2 quenched. However, also here it is difficult to establish the 
Pfaffian structure of the /c-point functions, and thus the number of integrals grows with k and Nf for 
the fe-point density correlation function (for recent progress see however |36j). 

It is an open problem how to include both Nf > and v > 1 in our skew-orthogonal polynomial 
method, in particular the latter, as the joint probability distribution of eigenvalues which we computed 
here for the general Nf and v case is no longer a simple Pfaffian of a single block matrix for v > 0. 
The reason for attempting this problem in the future is not only an independent confirmation for the 
aforementioned analytical results which are mainly for k = 1. The extension to arbitrary k is crucial 
if we want to get a hand on individual eigenvalue distributions at least in a perturbative expansion 
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(see e.g. [9]), apart from the conceptual aspect of integrability of our two-matrix model. Furthermore, 
in all works [20] [2Tj [26] [27] also the density of real eigenvalues of the complex non-Hermitian Wil- 
son Dirac operator was considered. The extension to this operator is another challenge to our method. 
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A Equivalence to the partition function of WchPT 

In this appendix we show the equivalence of our matrix model in the microscopic limit to WchPT in 
the epsilon regime, as it was displayed in eq. (j2.7|) . For non-degenerate parameters we have 



zi Nf) ~ / dAdBdWdVL JJ det 



(zf + m/)l n + A W + VL 

Wi + nt (z f -m f )t n+v + B 



g 2(l-a i ) 



* newt) 



Defining Zj = Zf±m,f we now express the determinants as integrals over Nf sets of complex Grassmann 
variables rjf , i = 1, . . . , n and ipl, b = 1, . . . , n + v, 

N f 

JTdet[...] = [ d 2 r]d 2 i)e v *^^ &ii+Aii)n i + ^^ w ^ , (A.2) 

j=i 

using summation conventions in all indices. A completion of the squares and integration of the 
Gaussian matrices yields 

= J dQ 1 dQ 2 dT det[Z+ - iV2(aQ 1 + T t )] n det[Z" - iV2(aQ 2 + T )]n+» e -Tv(TTl+Ql+Ql) 

(A.3) 

Here we have introduced two Hermitian matrices Q\p and one complex non-Hermitian matrix T 
of size Nf x Nf to do the Hubbard-Stratonovich transformation. Integrating out all Grassmann 
variables yields an exact expression containing the diagonal matrices with elements z^ (up to 
some constant prefactors that we dropped). If we parametrise T = UR with U unitary and R positive 
definite Hermitian the saddle point is taken at R ~ Ijsf.y/n when we employ the scaling of masses and 
a indicated in eq. (12. 9p . We obtain 



lim Z { u Nf) ~ f dQtdQ 2 ! ^det[[/]-ef r& ^ + ^- c/t ) +v/ ^ Tr ( ( 3 ll/+( 3^ t )^ riV W?+^) 

rwoo J J 

~ JdUdet[U} v exp ^Tr(M(U + + ^Tr(Z(U -U*)) - a 2 Tr(U 2 + ,(A.4) 

after integrating out the matrices Qi^ and performing a further rotation U — > ill, W — > —iU^. Eq. 
(|A.4[) is the partition function of WchPT in the epsilon regime for non-degenerate masses, hence we 
have completed our claimed proof. 
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A.l Alternative representation of Nf = 1 WchPT 



Finally we would like to make the equivalence between the Nf = 1 flavour partition function for the 
large- N limit of our even polynomials R^nix) explicit. The reason why we present this check in detail 
is that we require a slightly different form from the one given in [20J eq. (13), or in [21] eq. (151) (for 
z = there): 



-2a 2 



-2a 2 



d9 

— e^exp[mcos(6>) + iz sin(0) + 4a 2 sin 2 (6>)] 
2tt 



2vr 



dx 



7T 



e x exp[mcos(0) +izsm(9) — 4xasin(9)] , (A.5) 



where after using a trigonometric identity we have linearised the sine-squared term instead. We can 
now use the following identity [37] for generic matrices A, B 



L 



U{N f ) 



det[C/fexp 



Tt(AU + BU^) 



Nf-l 



2 N f (N f -l)/2 -Q j } 
3=1 



det[B]y /2 det^" 1 /^-! (//<)] 



det[A] 



(A.6) 

with /i? being the eigenvalues of the product matrix AB. Specifying to U = e %e G U(l) and identifying 
A = m + z + Aixd, B = m — [z + 4ixd) we arrive at 



zi N ^\m,z,d) 



-2a 2 



dx _ x 2 ( in — (I + 4ixd) ^ 
/% I m + z + 4ixa 



i//2 



j ^ \/ m 2 — (z + Aixa) 2 



(A.l) 



which agrees perfectly with eq. (|5.3p . see also eq. (99) in [27] . 



B Consistency checks of the partition function and factorisation 

In this appendix we perform a few consistency checks of the eigenvalue representation of the partition 
function eq. f)3. 15j) . For simplicity we will set m = and Nf = here. The antisymmetric weight eq. 
(I3.16P thus reduces to 



F(x) 



m=0 



2 e" 



"erf 



2V2C? 



(B.l) 



B.l The GUE limit a ->• 1 for v = 



On the level of the distribution eq. (|2.2p it is clear that the matrix W gets suppressed for a — > 1. 
After reducing to eigenvalues the same result should come out. Here we will look at the simplest case 
only, that is v = 0. We can insert the simplification eq. (|B.1|) into eq. (|3. 17|) at v = 0. Noting that in 
the limit a — > 1 the exponential factor in eq. (jB.ll) reduces to unity, we can use the following result of 
|28] (see end of section 3 in the first paper) for the limit of a Pfaffian of an error-function, 



lim Pf 

a— >1 



l<M<2n 



erf \ (di — dj 



'1-a 2 



A 2n ({d}) 



(B.2) 



This provides the second Vandermonde determinant in eq. ()3.17j) to constitute the GUE, and we 
arrive at 



lim Z 

a— >1 



(N f =0) l 
u=0 



\m=0 



N 



Y[ ddj d- f exp 



A N ({d}Y 



(B.3) 
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B.2 The chGUE limit a — > and factorisation 



Here we would like to perform a more detailed check. As it was shown in [21j in the large- N limit 
the spectral density of the Hermitian Wilson Dirac operator behaves like the density of a finite GUE- 
matrix of size vxv plus the microscopic chGUE density. We therefore expect that our general partition 
function will factorise in the limit a — > into a finite GUE times a second chGUE partition function 
of the remaining 2n eigenvalues. 

In order to show that we will first simplify Z^ further before using the de Bruijn integration 
theorem. Taking eq. (I3.13P as a starting point we can use the fact that the jpdf integrated over all 
variables is totally symmetric under exchange or relabeling of the d* @. We can split the variables 
into two sets, denoted by tj = \,...,v = dj—x,..., V i and sj=i J ...,2n = dj =v -\-i,„.,2n+v- If we do a Laplace 
expansion of the determinant into blocks of sizes 2n x 2n and v x u, we can split off a Vandermonde 
type determinant in the variables tj, after appropriately relabeling all other permutations: 



'3=1 



2n 



2 v 

J 



Y[dtje ^A u ({t}) / Y[dsje ^ II( s i ~ ^) A 2r i ({s}) 



(B.4) 



'i=i 



2a' 2 (l-a' 2 ) 



6=1 



e 2a z 



■ s 2n "1 

~~7T- 



. e 2a* 



1=1 

e 2a 2 . . . e 2a,2 



e 2a* 



e 2a 



2n "n 
1 - 



e 2a 



2m "1 
1 - 



1 h 
1 tu 



v-l 

1 

v-l 

V 



We can now perform the integrations over u±,. . . ,u n with the help of the standard de Bruijn inte- 
gral formula eq. (|3.14p . To make it applicable we can simply multiply the weight factors into the 
determinant, and we obtain the following result: 



2n 



7 " 



II dtj e~^A u ({t}f / Yl ds i Ife " ^) A 2n({4)Pfl<i J <2n [F( Si - Sj)].(B.5) 



1=1 



In eq. (|B.5p we can make further use of the symmetry under relabeling of the variables sj and the 
anti-symmetry of the Vandermonde determinant to write the Pfaffian as the following product: 



'i=i 



2n 



Yl dtj e ^A u {{t}f / Yl dsj e 3^ - U) A2n({s}) II F ^ 2k ~ S2k ~^- ( R6 ) 



i=l 



k=l 



Now consider the exponents of two consecutive variables, say si and s 2 , from the Gaussian and the 
exponential part of F: 



cxp 



s\ + s 2 2 
4a 2 



+ 



( S2 _ Sl )2 (1 _ a2) 



8a 2 



exp 



Qi + s 2 ) 

8a 2 



{si-s 2 f 



(B.7) 



In the limit a — > the first term on the right hand side will give a delta-function, whereas the second 
Gaussian term remains: 



lim V 8a 2 7r exp 



(si + s 2 ) 2 (si-s 2 ) 5 



6(8! + s 2 ) e" s ? /2 . 



8a 2 8 

In the same limit the error- function in eq. (jB.lj) becomes the sign function (see e.g. |28j), 



(Bi 



lim erf 

a->0 



(S2 ~ Si) 



T-a 2 
"So 2 " 



sign[s 2 - sij . 



(B.9) 



Using this symmetry our manipulations from now on will not be valid in general for correlation functions. 
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We are now ready to take the limit a — > of the partition function eq. (|B.6p , We first formally 
change variables rj = tj/a to isolate the a-dependence, and then take the limit only on the remaining 
a-dependence, that is the variables rj, 



lim 4°) 

a->0 



2 2\ 

a r,-, 



An({ S 2 }) 2 



r j J r S ■ 1 

U d rj e"- A,({r}) 2 / \[ d Sj SjSign[ Sj ] exp [J(sf 
>i=i ■ / -°° i =i i=i 

(B.10) 

We have used the fact that for S2fc — >• — S2fe-l we obtain A2 n ({s}) — >• Y\7=i A n ({s 2 }) 2 . The variables 
rj thus act as mass terms to the s-variables, but to leading order these masses can be neglected, leading 
to a complete decoupling of variables. Our final result is thus reading 



U dr 



■j exp 



A,({r}) 2 



2u+l 



exp 



A„({* 2 }) 2 , 



(B.ll) 



where the v would-be zero-eigenvalues at a = have decoupled into a GUE of size v, times the 
standard massless chGUE with exactly v zero eigenvalues. However, we should remember that due to 
rescaling the original variables tj are all of the order 0(a). 

It is clear how to this order of approximation the microscopic spectral density would look like: 
it is the superposition of a finite-^ semi-circle from the GUE and the Bessel spectral density with 
v zero-eigenvalues from the chGUE. This confirms the analysis of |21| on the level of the partition 
function, to that order. 



C Skew-orthogonal polynomials 

In this appendix we compute the skew-orthogonal polynomials given in eqs. fj4.6|) for v = 0. We 
first recall that the expectation value of a single characteristic polynomial (or equivalently the Nf = 
1 partition function) gives the even polynomials i?2j where we follow the argument of [38]. This 
expectation value is computed using the super symmetric method in the next subsection IC.ll Based 
on that the odd polynomials i?2j+i then follow as shown in lC.2l 

Starting from the quenched partition function eq. (I3.17P for v = 

/oo 2n 
H dSjW(Sj) A 2n ({ S }) Ph<i,j<2n[F( Si - Sj)] , (C.l) 

we can write the expectation value of a single determinant as follows: 

2n 



Y roc ^" 

(det[z - p 5 ]) 2n = — / Y[dsjw(sj)(z - Sj) A 2n ({s}) PhKij^niFisi - Sj)} 

j=l 

1 (2nY f°° N n 
= "ToT'Wi / H ds Msj) A 2n+1 ({ S },z) n^-^-i)- (C2) 

Zq ■ J -°° 3=1 1=1 



Here we have used the fact that both the Vandermonde determinant and the Pfaffian are antisym- 
metric in order to simplify their product (times the symmetric integral over all eigenvalues), cf. [29J 
Appendix 25. Furthermore we have included the determinant into a larger Vandermonde. Writing the 
Vandermonde determinant in terms of even and odd monic polynomials and ^2i+i respectively, 

A 2n+1 {{s},z) = det \{^ 2i { S j),^ 2i+1 {sj)} $2n{sj)] , (C3) 

t=0,...,n— l;j = l,...,2n+l L 
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and relabeling z = S2n+i, we have from the generalised de Bruijn formula (cf. [38 j eq. (4.9), or [34] 

(<Pi,<Pj) ■■■ fo(s2n+i) 



(det[ Z - $ 5 ]>2n 



z, 



(0) 



Pfi 



0<j,J<2n 



(C.4) 



^2n(s2n+l) 
-yo(s2n+l) • • ■ -¥ :> 2n(s2n+l) 

for the set of all monic polynomials <pi = {$0; ^3, ■ ■ ■ > ^2n}> with the scalar product defined 

in eq. (|4.9p . Choosing polynomials yjj = i?j to satisfy the skew-orthogonality eq. (|4.8p we obtained 
the desired relation 

R 2n (z) = (det[z - p 5 ]) 2n = (det[z + 5 }) 2n , (C.5) 

which obviously has parity. 

C.l Expectation value of a single characteristic polynomial 



We now compute the expectation value of a single characteristic polynomial eq. (jC.5ft using the 
super symmetric method for finite- N. As we have shown already this yields the subset of the even 
skew-orthogonal polynomials R 2n {x) from which the odd polynomials are constructed. At the same 
time this also gives the one flavour partition function, which we can compare to the known result from 
Wilson chiral perturbation theory [20] in the large- N limit as a consistency check. 

The calculation we present here for v > is slightly more general than needed in eq. (14.61) . Using 
the definitions (|2.ip and (|2.1ip for the Hermitian matrices A and B of sizes nxn and (n + v) x (n + v) 
as well as matrices W, Q, of sizes n x (n + v) with weights eqs. ([2. 2D we have for N = 2n + v 



(det[z + 5 ]> 



1 



A" 



2, 



(0) 



Xe 4a 



dAdBdWdfl det 



^Tr^+B^nQt) 



Wt + nt («-m)l n+ „ + B 



(C.6) 



Expressing the determinant through integrals over complex Grassmann variables rji, i = 1, . . . ,n and 
t(3b,b = l,...,n + u, 

det[...] = J # T j#^ e Vi((*^) s v+ A v)vjW>Uwh+^ ^ ^ c _ 7 ^ 

we obtain the following result after performing the Gaussian integrals over the matrices A, B, W, f2: 
(det [2 + 5 ])tv ~ J d 2 r]d 2 ifj e ( z+m HvM^-rn)ij;A-a 2 {(v*v l ) 2 H^i>b) 2 )-2^ b i> b v*v t 

d 2 U 



d 2 r]d 2 ip 



d,t 



TT 



~* |u| exp[(z + m - 2iat)rj*r]j + (z + m - 2iat)iplip b ] 



x exp[-i\/2(l - a 2 )(u7]*rij + uipl^b)} 



dt 



x z — m 



d 2 u 



TT 



-* 2 -|«| 2 



+ m + 2iat + i\/2(l - ° 2 ! 



+ 2iat + 1^2(1 - a 2 )u* 



n+v 



(Ci 



In the second step we have linearised upon using one real and one complex variable in our Hubbard- 
Stratonovich transformation. As a last step we use the following integral representation for generalised 
Laguerre polynomials L^(x) [39] 

— e"!"! 2 (iu + y) n (iu* + \) n+v = (-) n A"4" ) ( A 7) , (C9) 
/c n 
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to obtain 



(det[z + 5 ]> 



N 



n!(2(l -a 2 )) n+ 2 



00 ~~-*( Z + 2iat - m YrM 



[z + 2iatf 



m 



2 1- a 2 



.(C.10) 



Here we have put the correct normalisation coefficient that can be read off from (detjz + 05])at = 
z N + 0{z N ~ 1 ) for large z. For v = this reduces to the result given in eq. (|4.6p . It can be easily seen 
that in the microscopic limit eq. (15. ip using eq. (15. 2ft this reduces to the Nf = 1 flavour partition 
function eq. ()A.7j) 



C.2 The odd polynomials R 



2j+l 



X) 



In this subsection we construct the odd set of skew orthogonal polynomials directly from the even 
ones by differentiation. Instead we could have computed the R2j+i{x) again from the supersymmetric 
method as in the previous subsection, using the relation [38] R2 n +i(x) = (det[x — Ps\(x + Tr05))2 n - 
Based on our result above eq. (IC.lOj) for v = K% n {x) = (det[x — ]p^\)2n we define 



R2j+i(x) 
w(x) 



-2a 2 w(x)~ l — 
ax L 

exp[-x 2 /(4a 2 )] . 



w(x)R2j(x) 



(C.ll) 
(C.12) 



In the following we will show by direct computation that these satisfy the skew-orthogonality relations 
eq. (14. 8p with respect to the scalar product 



(h,g) 
F(x) 



dxdy w{x)w{y)F {x — y)h{y)g{x) 



(C.13) 



e Sa^ 



xy/i - a 2 
erf ( ■ h 
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Y^dHl 



+ erf 



x 



W ^/2a 2 {\ - a 2 ) 



Because the polynomials obviously have parity, Rj(—x) = (—) J Rj(x), and because the weight is 
antisymmetric, F(—x) = —F(x), it holds 



(R2j,R2i) — — (R 2 j + i, R21+1) 



(C.14) 



Furthermore, the even polynomials are skew-orthogonal to all polynomials of lower degree, as can be 
seen from eq. (|C4[) , 

(R2j,R2i+i) = , for j>l. (C.15) 
It remains to investigate the remaining cases. For j < I we have 



(i?2j> R%H 



2a J dxdyw(y)F(x - y)R 2 j(y)— w(x)R 2 i(x) 

f dF(x — y) 
-2a j dxdyw(x)w(y) — R 2 j(y)R2i(x) 



dy 



d 



2a j dxdyw(x)F(x - y)R 2 i(x)— w{y)R2j{y) 



After a first integration by parts with respect to x we have used dF ^ x y ^ 



(C.16) 



dx 



integration by parts with respect to y yields zero after using eq. (|C.15|h 



■ dF ^ y , and a second 
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It remains to compute the only nonvanishing product for j = I: 

/dF(x — y) 
dxdyw(x)w{y) — R 2j (y)R 2 i (») 



1 - a 2 

dxdyw(x)w(y) — - — (x - y)F(x - y)R 2j (y)R 2l (x) 



jQ l Qp'(\ — — 771,2 f ( rn{x — y) —m(x — y)\ 

+ \l -e taHi-o*) / dxdyw{x)w{y) ( e 2^ + e 2^ j R 2j (y)R 2 i(x) 

a 2 (l-a 2 ) J dxdy[^w(y)-w(x)^)F(x-y)R 2 ,(x)R 2j (y) 



2o 2 (l - a 2 ) 



„2 



7T 



+2\ -e 2 ( 1 - a ") / dxw(x - m)R 2j (x) / dyw(y + m)R 2j (y) . (C.17) 



In the first step we simply integrate by parts. We then observe that the first line in the last equation 
is proportional to (R 2 j, R2j+i) because when the derivatives also act on the polynomials R 2 j(x) and 
R2j{y) this gives zero as the polynomials R 2 j are skew-orthogonal to all lower order polynomials. 
Putting this term on the left hand side and using parity we obtain 



a 2 {R 2j ,R 2j+1 ) = 2^ 2q2 ^ a % (J dxw(x - m)R 2j {x)\ . (C.18) 

In a last step we have to evaluate the integral over the weighted even polynomials. Using eq. (jC.lOj) 
we have 

f r , \ t~> f \ j!(2(l-a 2 )V' f°° , j*-™.) 2 [°° ,0 f(x + 2iat) 2 - m 2 \ 
jdx W (x-m)R 2] (x) = J jxe ^J_Jt^L ^- — J^- j 

j!(2(l-a 2 )V 1 f°° , , -l+£ T fx 2 + 2isx- s 2 + 2m(x + is)\ 

= ; — . . ^_ / dxds e 4» 2 Lj : ^ 

C-^VtF 2a j \ 2(1 -a 2 ) J 

j\ (2(1 - a 2 )y 



2avr . (C.19) 

In the last but one line we can change to complex coordinates z = x + is and to an integration over the 
full complex plane. Because of the argument of L n (z(z + 2m)/2(l — a 2 )) the angular integration will 
give zero for all powers of the argument z, except for the constant term z° which is unity, L n (0) = 1. 
Collecting all factors we thus obtain 

Tj = (R 2j ,R 2j+ i) = 8^a2 2j (j!) 2 (l-a 2 ) 2j+5 e . (C.20) 
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